Power functional theory for Brownian dynamics 
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For classical many-body systems subject to over-damped Langevin dynamics, we propose a dy- 
namical functional consisting of reversible free energy changes and irreversible power dissipation. 
Minimization with respect to the microscopic one-body current yields a closed equation of motion. 
Dynamical density functional theory is obtained when using the ideal gas power functional. Ap- 
proximations to the excess (over ideal) power functional yield numerically tractable equations of 
motion beyond the adiabatic approximation, opening the door to the systematic study of systems 
far from equilibrium. 
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The dynamics of soft matter is a topic of considerable 
current interest, both from an experimental and theoreti- 
cal perspective [U [2] . Among the various theories aiming 
to provide a microscopic description of complex trans- 
port and pattern formation phenomena, the dynamical 
density functional theory (DDFT) [3J 0] has emerged as 
a prime candidate. Much of the appeal of DDFT arises 
from its close connection to equilibrium density func- 
tional theory (DFT), which is a virtually unchallenged 
framework for the study of equilibrium properties in in- 
homogeneous liquids [3J. 

The DDFT was originally suggested on the basis of 
phenomenological reasoning [3J, and has since been red- 
erived from the Langevin equation 0], using projection 
operators [5] and via coarse graining the many-body 
Smoluchowski equation j6j. Despite the deep insight 
provided by these alternative derivations of the DDFT, 
the final equation of motion for the density remains un- 
changed and no guidance is provided for making a sys- 
tematic improvement that goes beyond the standard for- 
mulation. A striking and undesirable feature of DDFT 
is the clear asymmetry between the treatment of spatial 
and temporal degrees of freedom. The intricate nonlo- 
cal spatial structure is in stark contrast to the simple 
time-local treatment of the dynamics and suggests much 
potential for improvement, particularly in view of the 
sophistication of memory function approaches such as 
mode-coupling theory [7]. 

At the heart of DDFT lies the 'adiabatic approxima- 
tion', in which the nonequilibrium pair correlations of 
the real system are approximated by those of a fictitious 
equilibrium, whose density distribution is given by the 
instantaneous density of the nonequilibrium system [5]. 
Although reasonable for describing relaxational dynamics 
close to equilibrium, this assumption fails when applied 
to even the simplest driven systems (e.g. shear flow [3]), 
thus putting out of reach many relevant and fundamental 
non-equilibrium problems. 

In this Letter, we overcome the above limitations by 
developing a variational theory based on the current 



rather than the density as the fundamental variable. 
Our framework transcends the adiabatic approximation 
and provides a physically intuitive method to go beyond 
DDFT. 

The central variable in DDFT is the time-dependent 
one-body density, 

p(r,t) = (^2S(r-r i (t))y (1) 

where (•) indicates an average over trajectories, Vi(t) is 
the position of particle i at time t, S(-) is the Dirac distri- 
bution, and the sum is over all N particles in the system. 
However, the motion of the system is better characterized 
by the time- and space-resolved one-body current 

J(r,t) = ^v i (t)«(r-r < (t))^ ) (2) 

where Vj(t) = dri(i)/dt is the velocity of particle i. It 
is evident that the vector current ([2| provides a more 
appropriate starting point than the density for a the- 
ory aiming to describe flow in a complex liquid; there 
are many relevant flows for which the density is trivial, 
p{v,t) = const, but the current is not. Once the current 
is known, the density can be calculated from the conti- 
nuity equation, 

p(r,t) = p(r,t )- f dt'V ■ J(r,t'), (3) 

J t 

which expresses the local conservation of the particle 
number. Prior to the initial time t we take the system 
to have been in equilibrium. 

Starting from a particle-based description, we write the 
(over-damped) Langevin equation as 

7 v 2 - Fj = 0, (4) 

where 7 is a friction constant. The viscous force is bal- 
anced by the total force 

F, = -ViU(r N )- V t V cxt (v t ,t)+X(v t ,t) + xM, (5) 
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consisting of four contributions: (i) Forces arising from 
the total interaction potential, U(r N ). (ii) Conservative 
forces generated by an external potential, V ex t( r >t). (iii) 
Non- conservative forces X(r, t). (iv) Stochastic, ther- 
mal forces, Xi(t), which are auto-correlated according 
to (Xi(t)Xj(f )) = 2k B T-/6 tj 5{t - t')l, where k B is the 
Boltzmann constant and T is temperature. 

In the spirit of the Rayleigh dissipation function [T0] , 
but including both dissipative and conservative forces, 
we define a generating function 



N V N 



>-£( 
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(6) 



obtained by integrating the left hand side of Q for- 
mally with respect to Vj and summing over all particles. 
Clearly, the Langevin equation Q is then given by 



dR(r f 



0. 



(7) 



where the derivative is taken at fixed particle positions, 
r N = const; the physically realized velocities v N min- 
imize R. When averaged and integrated over time the 
function R plays the role of a generating functional for 
the density ([I]) and current ^ via 

p(r,t) = - . S —— I dt(R), (8) 



J(r,i) 



5V ext (r,t) 
8 



dt(R), 



(9) 



where V ex t(r,t) = dV exi (r,t)/dt. We can hence obtain 
the intrinsic part of the generator J dt(R) as a functional 
of p(r,t) and J(r, t) via Legendre transform 



W[p,3] 



dt(R) 



(10) 



+ Jdtjdr (X(r,t)-J(r,t) + V«t(r,t)p(r,t)), 

which we split into a sum of dissipative and reversible 
contributions, 



W[p,3} = J dt (P[ P) J}+F{ P } 



(11) 



Here the integrand consists of a dissipated power func- 
tional, P[p, J], accounting for irreversible energy loss due 
to the friction, and the total time derivative (indicated by 
the dot) of the intrinsic Helmholtz free energy functional, 
F[p], responsible for reversible free energy changes [3]. 
In general, the dissipation functional P[p, J] will be non- 
local in space and time, depending on the history of the 
fields p(r,t) and J(r, t) prior to time t. Employing the 
functional chain rule, F[p] = fdr p(r,t)6F[p]/8p(r,t), 
and the continuity equation ph, followed by a partial 
integration in space leads to 



W[p,3] =j dt (p[p,3] + J drJ(r,t) 



8p(v,t) 



(12) 



This intrinsic contribution can be combined with the ex- 
ternal work in order to define a total dissipation func- 
tional 

n\p,3]=W\p,3] (13) 
- J dt J dv (X(r, t) - W oxt (r, *)) ■ J(r, t). 

Recognizing Q, an equation for the time evolution of 
the system is obtained [IT] by 



sn\p,3\ 

<yj(r,i) 



(14) 



at p(r, t) = const. The result is our fundamental equation 
of motion 



(15) 



where the temporal variable in |llj) has been relabelled 
if , and we used the shorthand notation p' = p(r,t') and 
J' = J(r,t'); furthermore SP[p' ,3'}/S3(r,t) = for t > if 
due to the fact that P[p',3'] depends only on the value 
of its arguments at times prior to if, cf. the discussion of 



this causal structure below ( 11 ) 



The term on the left hand side of ( 15 ) represents a fric 



tion force, balanced by the terms on the right hand side 
arising from inhomogeneities in the local intrinsic chemi- 
cal potential, 8F[p]/5p(r, t), and external forces. In equi- 
librium the left hand side of ( 15 ) vanishes and the vari- 
ational prescription reduces to 8Cl[p]/5p(r) — 0, where 
tt[p] is the grand potential functional [3], thus recovering 
equilibrium DFT as a special case. In nonequilibrium, 
the challenge is to find explicit forms for the dissipation 
functional, P[p, J], to obtain a closed equation of motion. 
To gain some intuition we consider three special cases: 

(i) Frictionless motion. The simplest case is obtained 
by neglecting nonconservative external forces and setting 



P[p,3] = 0. Equation Q then yields SF[p]/5p(r,t) = 
P — V cyL t(r, t), which is identical to the Euler-Lagrange 
equation of equilibrium DFT, where the chemical poten- 
tial, p 7 keeps the particle number constant. In this limit 
the theory mimics "ultrafast" dynamics, where the den- 
sity field instantaneously follows changes in the external 
potential. 

(ii) Ideal gas. For a system of nonintcracting parti- 
cles the Helmholtz free energy is known exactly, -Fid [p] = 
k B T j drp(r, t)[ln(p(r, t)A 3 ) - 1], where A is the thermal 
wavelength, and the dissipation functional is given by 



Pid[p,3] 



dr 



7J(r,t) 2 
2p(r,t) 



(16) 



Functional differentiation of this expression with respect 
to the current generates the (mean) friction force, such 
that J(r, t)-8P[p, 3]/S3(r, t) = 3(r,t)"yv(r,t) represents 
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the dissipated power density resulting from the average 
motion of the system; here v(r, £) = J(r,t)/p(r, t) is the 
average velocity. The exact expression for the ideal cur- 
rent, 7 J id (r, t)/p(r, t) = -k B TV In p(r, t) - W cx t(r, t) + 



X(r, i) follows from substitution of (16) into (15). 

(in) Dynamical density junctional theory. The intrin- 
sic free energy functional of an interacting system can 
be written as F[p] = F;d[p] + F exc [p], where the excess 
contribution accounts for interparticle interactions. If we 
retain the full free energy, but assume that the dissipation 
is given by (16), then (15) yields J(r,i) = Jddft^ 1 ") i), 
where 



7 J 



DDFT 



M) 



P(r,t) 



, SF\p] 
Sp(r,t) 



W ext (r,t) + X(r,<) (17) 



is precisely the current of DDFT [31 0]. We thus gain 
new insight into the standard theory, namely that the 
adiabatic approximation is equivalent to assuming a triv- 
ial, noninteracting form for the dissipation functional. 
This observation suggests that superior theories can be 
obtained by developing approximations to P[p, 3] which 
recognize the existence of interparticle interactions. 

We can now go beyond DDFT by decomposing the 
dissipation power functional into two contributions, 



P[p,3]=P id [p,3}+P cxc [p,3}, 



(18) 



where Pjd [p, 3] is given by ( 16 ) and P cxc [p, 3} accounts for 
dissipation arising from interparticle interactions. Mi- 
croscopically, the excess dissipation is caused by corre- 
lated fluctuations of the particle velocities Vi(t) around 



the mean velocity v(r, ; (i), t). Using ( 18 1 the fundamental 
equation of motion (fl5|) becomes 



7-JM) 



dt 



^Pexcfp'.J'] 



<JJ(r,t) 



6F[p] 



W«t(r,t)+X(r,t). (19) 



5p(r,t) 

Corrections to DDFT thus come from the second term 
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FIG. 1. Free relaxation of a density peak (slab geometry) 
from an initial state (red dashed curve) for which the density 
is confined to a parabolic potential. The solid green curves 
show the density and the blue dash-dotted curves show the 
spatially resolved dissipation (scaled down by a factor 20) for 
three different approximations to the excess dissipation (at 
dimensionless times t* = tksT '/yd 2 — 0.2 to 1.4 in steps of 



0.2). a) Using DDFT, b) using p4| , c) using (26]). Panel d) 
shows how the total dissipated power decays as a function of 
time. 



with ( 19 ) yields the identification 
,6P„ 



dt'- 



83(r,t) 

SF CKC [p] 
5p(r,t) 



(21) 



o^lr r' t) 
dr' P [ ' J V4>(r,rO. 



P(r,t) 

If we split the full nonequilibrium two-body density into 
an instantaneous-equilibrium and an irreducible part, 
/9 (2) (r,r',t) = pcq (r, r',£)+/5; rr j (r, r',i), then we can iden- 
tify the following two identities 



dt' 



SF exc [p] 
8p(r,t) 

SPcxcIp',3'} 
JJ(r,t) 



dt 



= / dr 



V^(r,r'), (22) 



(23) 



Equation ( 22 ) is an exact equilibrium relation [3] , 



on the left hand side. It is instructive to compare (fl9|) whereas (|23j) expresses the fact that the excess dissipa- 



with the formally exact result obtained by integrating 
the many-body Smoluchowski equation over N — 1 par- 
ticle coordinates [5]. Restricting ourselves for simplic- 
ity to systems that interact via pair potentials, U — 
^2i^tj 0( r i: r j)i an exact equation of motion is obtained, 



7J(r,t) 

p(T,t) 



fc B TVln( / 9(r,t)A 3 )-VKKt(r,t) 
,P (2) (r,r>,t) 



X(r,t)- J dr' 



P(r,t) 



V^(r,r'), (20) 



tion is intimately connected with beyond-adiabatic, ir- 
reducible two-point correlations. This is fully consistent 
with the interpretation that the dynamic functional is a 
generator for many-body induced friction forces in the 
system. 

Any practical application of the framework, i.e. imple- 
mentation of (19) along with ([3]), requires one to choose 
an approximation for P exc [p, 3] for the given physical sys- 
tem (as characterized by U). Here we discuss several 
model approximations to P eX c, the simplest of which is 
given by the space- and time-local expression 



where p( 2 '(r,r' ,t) is the exact (but unknown) out-of- 
equilibrium equal-time two-body density. Comparison 



^°c C) [/M]= / dr (7(p)-7) 



JM) 2 
2p(r,ty 



(24) 
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where j(p) is a density-dependent friction factor which 
can be approximated using either a virial expansion or, 
in the case of hard spheres, by employing closed form 
expressions [IS]. Using (24) slows the dynamics in re- 



gions of high density relative to DDFT [T3] , but will still 
fail to describe collision induced dissipation arising from 
shearing motions in the fluid. To capture this requires a 
nonlocal treatment in space and in time which recognizes 
the finite range of the colloidal interactions. A suitably 
general form is given by 



p(nloc; 
exc 



[p, J \=J dr J dr'Jdt' J(r, t) ■ K(r, t; r', t') ■ J(r', t') 



(25) 



where the convolution kernel K is a second-rank tensor 
describing current-current scattering and is in general a 
functional of p. The possibility of incorporating temporal 
nonlocality represents a powerful feature of the present 
approach. A simple specific form for the scattering ker- 
nel which captures this is K = 7m(£ — £')<5(r — r')//9(r, i), 
where m(t — i') is a memory function (with units of in- 
verse time). As a result we obtain a spatially local excess 
dissipation functional with temporal memory given by 



Pt cla) [ P J}= f dr f 

J Jtn 



yj(r,*Mfro m(f _ 



(26) 



For general interactions and time-dependent external 
fields the memory need not be time translationally in- 



variant, as assumed here. The equation of motion ( 19 ) 
thus becomes 



(J(r,t)+ / dt'3{r,t')m{t-t') 



, SF[ P ] 
6p(r,t) 



WextM)+X(r,i). (27) 



We have performed numerical calculations using the 
(three-dimensional) hard sphere system in a slab geom- 
etry using approximations (24) and (26). The (spatially 



non-local) Rosenfeld functional [14] was used to approx- 
imate -F e xc[/>]- The system is initially confined along the 
y-axis by a parabolic potential, V ex t(r,t < 0) = by 2 /2, 
where b is a constant measuring the strength of the trap- 
ping potential. We then switch off the confinement at 
t = to = 0, such that Kxt(r, t > 0) = 0, and monitor the 
time evolution of the density. In Fig. [I] we show the den- 
sity calculated using three approximations: DDFT, the 
density dependent friction coefficient ( p4| (where we have 
employed the expression for 7(p) from [12]) and, finally, 
the temporally nonlocal approximation ( 26 1 . When im- 



plementing ( 26 ) we have approximated the memory func- 
tion by the simple form m(t—t') — (a/r) exp(— (t— t')/r), 
where a is a dimensionless parameter and r is a relax- 
ation time. 



The density profiles in Fig. [I] (generated using the pa- 
rameter values a = 0.8,6 = l0k B T,r = 0.05jd 2 /k B T, 
where the particle diameter, d, sets the length-scale) 
show that both (24) and (26) slow the dynamics. The 



magnitude of the retardation achieved using ( 26 ) depends 



upon the choice of parameters a and t determining the 
memory. While ( 26 ) generates density profiles with a 
similar functional form to those from DDFT, Eq. (24) 



leads to a density more sharply peaked at the origin. The 
spatially resolved dissipation (blue dash-dotted curves) 
confirms the intuition that power is dissipated mostly 
in regions of high density gradient and the spatial inte- 
gral of this quantity, J dr J(r, t) ■ SP[p, J]/(5J(r, t), decays 
towards zero as the system approaches equilibrium. It 
is well known that relaxation rates predicted by stan- 
dard DDFT are significantly faster than those found in 
simulation [H [TJ] . On the basis of our findings we con- 
clude that this failing can be remedied by the incorpo- 
ration of temporal nonlocality in the excess dissipation 
functional, leading to memory effects in the equation of 
motion for J(r, t). If an exponential memory is employed 
then the computational demands are comparable to those 
of DDFT. 

In summary, we have shown that collective Brownian 



dynamics can be formulated as a variational theory ( 14 ) 



based on the dissipative power as a functional of the one- 
body density and the one-body current. The generaliza- 
tion to mixtures of different species is straightforward: 
the ideal contribution becomes '^2 i Pid[pi, Ji], where i 
enumerates the different species, whereas P e ^ c [{Pi}, {Ji}] 
generates dynamical coupling between particles of differ- 
ent species [IB] . 

The appeal of our approach stems from its utility and 
and ease with which it can be implemented, which con- 
trast strongly with the time-dependent DFT formalism 
by Chan and Finken [IB]. As far as we are aware, the 
approach of [16] has never been applied to any model 
system, possibly because it is built around an action func- 
tional, which could not (yet) been approximated in any 
systematic or physically intuitive way. This is not the 
case with the excess dissipation functional identified in 
the present work, for which even simple expressions, such 
as ( 26 ) , transcend the adiabatic approximation. While 



DDFT remains an active field of research (e.g. hydrody- 
namic interactions |17] and arbitrary particle shapes |18] 
have very recently been addressed), it is important to ap- 
preciate that all extensions and modifications proposed 
since the original presentation of DDFT [3] have been 
firmly under the constraint of adiabaticity. The power 
functional approach is free of this restriction. 

Finally, the relationships ([8| and (|9| of the one-body 
fields to their generating functional imply that two- and 
higher-body dynamic correlation functions can be gen- 
erated from further functional differentiation, putting a 
nonequilibrium Ornstein-Zernike relation within reach. 
Work along these lines, as well as applications to driven 



5 



lattice models [TS], is currently in progress. 
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